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ABSTRACT 

We present the largest sample to date of giant molecular clouds (GMCs) in a substantial spiral galaxy other 
than the Milky Way. We map the distribution of molecular gas with high resolution and image fidelity within 
the central 5 kpc of the spiral galaxy NGC 6946 in the 12 CO (J = l-0) transition. By combining observations 
from the Nobeyama Radio Observatory 45-meter single dish telescope and the Combined Array for Research 
in Millimeter Astronomy (CARMA) interferometer, we are able to obtain high image fidelity and accurate 
measurements of Leo compared with previous purely interferometric studies. We resolve individual giant 
molecular clouds (GMCs), measure their luminosities and virial masses, and derive Xco - the conversion 
factor from CO measurements to H2 masses - within individual clouds. On average, we find that Xco — 1.2 x 
10 20 cirT 2 / (K km s _1 ), which is consistent within our uncertainties with previously derived Galactic values as 
well as the value we derive for Galactic GMCs above our mass sensitivity limit. The properties of our GMCs 
are largely consistent with the trends observed for molecular clouds detected in the Milky Way disk, with the 
exception of six clouds detected within ~400 pc of the center of NGC 6946, which exhibit larger velocity 
dispersions for a given size and luminosity, as has also been observed at the Galactic center. 
Subject headings: ISM: molecules - galaxies: ISM - galaxies: individual (NGC 6946) 



1, INTRODUCTION 

Studies of extragalactic giant molecular clouds (GMCs) 
shape our understanding of spiral galaxies on both small and 
large scales. Within GMCs, on scales smaller than these 
massive bound gas structures, individual gas cores inside the 
clouds condense and form new stars. On scales larger than 
GMCs, global galaxy dynamics govern the motions of these 
massive gas clouds, where we study the physical environ- 
ments which create and destroy them. Much of what is known 
about the evolution of the interstellar medium in spiral galax- 
ies has been learned by analyzing the propert i es of these 
large molecular clouds of gas (i.e.. [Wa da 2008; [Koda et al.l 
2009; iTasker & Tanl 120091 iDobbs. Burkert. & Pringld 1201 it 
Egu sa, Koda. & Scovilldl201 ll) . 

In order to study the gas of GMCs, which is primarily 
composed of molecular hydrogen, other molecular tracers - 
the most common of these being the lower rotational tran- 
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sitions of CO - are typically observed since the overall gas 
temperatures are too low to directly excite H2 line emission. 
Accurate determinations of the masses of GMCs are thus 
dependent upon an accurate relation between CO flux and 
H2 mass, known in the literature as the Xco factor. Val- 
ues of this conversion factor have be en derived for molec- 
ular clouds within our own Galaxy (ISolomon et al.1 119871: 
IScoville et al.lll987t iDame. Hartmann. & Thaddeusll2001l) a s 
well as compiled for nearby galaxies (Bo latto et al.l 2008): 
to within a factor of two, the typic al value seems to ho ver 
aro und 2-3 x 10 2Q cm" 2 / (Kkm s' 1 ) (ISolomon et al.ll987l but 
seelHever et al.ll2009t[Scoville et al.lll987l;IWilson & Scovilld 
19901; IDame. Hartmann. & Thaddeusl 120011; lOka et al.l l2001t 
Blitz et al.H2007t iBolatto et al.ll2008l). with higher values ob- 
served in the SMC ([Israelii 19971: iLerov et al.ll2007 l lBlitz et all 
l2rj0^lBolatto^FaIl l2008). 

The conversion factor, derived from Galactic clouds and re- 
solved extragalactic GMCs in various environments, is used to 
derive molecular gas masses for unresolved GMCs in galaxies 
out to high redshift. Deriving this value in nearby GMCs is 
therefore crucial to our understanding of the amount of molec- 
ular gas, molecular fraction, and star formation efficiency (star 
formation rate per unit gas mass) in galaxies as a function of 
redshift and environment. 

For galaxies other than the Milky Way, GMC stud- 
ies are typically resolution -limited (see the discussion in 
iRosolowskv & L erov 2006). In order to identify individ- 
ual molecular clouds, only nearby galaxies can be stud- 
ied, which has kept the available sample of extragalac- 
tic clouds relatively small. For this reason, investigations 
into individual extragalactic GMC properties, and varia- 
tions of these with GMC environments, have been mostly 
limited to clouds inhabiting Local Group galaxies (e.g. , 
[ Wilson & Scovilldll990tlArimoto. Sofue. & Tsui imot oll 19961; 
l lsraell 119971; iBoselli. Lequeux. & Gavazzil 120021; iBlitz et al. | 
l2007t iBoiatto et al.1 l2008t iHughes et al.l 120 10t iBigiel et all 
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l20l0l:lLerov etalEoTlh . 

Even when they are resolvable, the determination of GMC 
boundaries is not straightforward. In addition to being lim- 
ited by instrumental resolution, the appropriate spatial separa- 
tion of two (or more) adjacent clouds is not always obvious in 
crowded or highl y blended regions. Var i ous studies of Galac- 
tic GMCs (e.g., iSolomon et all [19871 iScoville et alJ fl987l) 
have solved this issue by identifying clouds "by hand" above 
a specific temperature contour (which is dependent upon the 
sensitivity of the instrument). Mo re recently, automated algo- 
rithms have been developed (e.g.. [Williams, de Geus. & Blitz 
Il994t iRosolowsky & Leroyl 120061) in order to identify indi- 
vidual clouds; these commonly used codes differ in their 
clump identification philosophy, particularly with regard to 
the amount of flux included in identified clumps, as discussed 
at length by IRosolowsky & L eroy (200g). 

Before even identifying individual GMCs, many studies of 
extragalactic clouds suffer from the more fundamental prob- 
lem of missing zero spacing information. Observing the gas 
distribution of a galaxy with only an interferometer neglects 
extended flux on the largest scales, where it is "resolved out" 
by the beam of an interferometric array. This effect can be 
quite significant (Ko da et al.ll201 ll) . 

1.1. CARMA &■ NR045 CO Survey of Nearby Galaxies 

To investigate GMC evolution in galactic disks and re- 
solve the physics which controls the star formation rate within 
GMCs, we will present resolved observations of GMCs across 
the disks of nearby galaxies in a CO Survey of Nearby Galax- 
ies based on observations using the CARMA interferometer 
and the Nobeyama Radio Observatory 45-meter (NR045) sin- 
gle dish telescopes. To date, 13 galactic disks have been ob- 
served with both CARMA and NR045 as part of the survey, 
and 4 additional galaxies have been observed with the single 
dish. The design of the survey will be discussed in Koda et 
al. (in prep). In our survey, we will resolve individual GMCs 
in a significant sample of nearby spiral galaxies with a variety 
of morphologies in order to study the evolution of molecular 
clouds and star formation. 

In this paper, we combine observations from CARMA 
and NR045 in order to achieve high resolution, as well as 
extremely high image fidelity, and resolve individual ex- 
tragalactic GMCs in NGC 6946 to highlight the results 
made possible by our survey. In Section II we describe 
the observations, and we present the results in Section III 
as well as a discussion of the boundaries of individual 
GMCs. We utilize the well known CLUMPFIND algorithm 
d Williams, de Geus. & Blitzj 19941) to derive sizes and velocity 
dispersions for individual GMCs. In Section IV, we discuss 
the properties of our detected GMC sample and determine the 
conversion factor Xco within individual clouds. We summa- 
rize our conclusions in Section V. 

2. OBSERVATIONS 
2.1. Nobeyama Radio Observatory 

We observe NGC 6946 in the 12 CO ( 1 -OVtransition with the 
Nobeyama 45-meter single dish telescopq3> using the Beam 
Array Receiver System (BEARS) instrument. Observations 
were performed during the early months of 2008, 2009, and 

12 Nobeyama Radio Observatory is a branch of the National Astronomical 
Observatory of Japan, National Institutes of Natural Sciences. 



2010 (throughout our three year observing program at NRO) 
as part of our CO survey of nearby galaxies. BEARS is a 
multi-beam receiver with 25 beams, which are aligned in a 
5x5 orientation. The FWHM of the 45-m dish is 15" at 
1 15 GHz (19.7" after regridding), and we observe with chan- 
nel increments of 500 kHz and Hanning smooth for a velocity 
resolution of 2.54 km s _1 . After dropping the edge channels, 
we use a bandwidth of 265 MHz (690 km s ). For most 
of the scans, the system temperature in the double side band 
(DSB) of BEARS ranges from 300-400 K; scans with T^ 
much higher than 400 K (due to being observed during day- 
light hours) are heavily down-weighted and as a result do not 
contribute much to the final images. The ratio of the upper- 
to-lower side band (known as the scale factor) was confirmed 
each year to be within a few percent by observing Galactic 
CO sources (e.g., AGB stars) with BEARS and the single 
side band (SSB) S100 receiver and making subsequent cor- 
rections. We use the S40 receiver and Galactic masers for 
pointing. The pointing was checked roughly every two hours 
during the observations and was accurate to 2-3". We con- 
vert Ta* to T m h assuming that the main beam efficiency of the 
telescope is 0.4 (i.e., T m j = Ta- / 0.4). 

Using on-the-fly (OTF) mapping, NGC 6946 was scanned 
in the RA and Dec directions, and positions external to the 
galaxy (OFF positions) were observed between scans. Ex- 
trapolating between OFF scans on opposite sides of the galaxy 
greatly reduced non-linearities in the spectral baselines. The 
duration of each scan (ON + OFF) was ~1 minute, and the 
entire galaxy was mapped in ~40 minutes; a total of 36 us- 
able maps were taken for a total of 24 hours of observation 
time (including ON, OFF, and slew time). The scans were 
separated by 5", resulting in oversampling by a factor of 3 
compared to the 15" FWHM of the beam, which is neces- 
sary in order to achieve Nyquist sampling (5.96") of AcofD 
= 11.92" (where Xqo is the observed wavelength, 2.6 mm, 
and D is the antenna diameter, 45 m). The data reduction 
and sky subtraction, performed by interpolating between the 
OFF scans for each OTF (ON) scan, were completed using 
the NOSTAR package developed at the Nobeyama Radio Ob- 
servatory. Spatial baseline-subtracted maps were made sepa- 
rately from the scans in the RA and Dec directions in order to 
minimize systematic errors in the scan directions, and these 
were subsequently co-added. The rms noise of these single 
dish observations is 0.13 K (0.57 Jy beam -1 ). 

2.2. CARMA 

To complement the single dish observations, NGC 6946 
was also observed in the CO (1-0) transition in April 2009 
with the C and D configurations of CARMA. CARMA is a 
15-element interferometer which combines six 10-meter an- 
tennae (originally the Owens Valley Radio Observatory, or 
OVRO) with nine 6-meter antennae (formerly the Berkeley- 
Illinois-Maryland Association, or BIMA) to achieve superior 
uv-coverage compared to either of its predecessors. The ob- 
servations of NGC 6946 were performed using three dual side 
bands, each with 63 channels, for a total bandwidth of ~100 
MHz (after removing six edge channels per sideband) and a 
channel width of 2.54 km s . After a total of ~21 hours on 
source (including calibrators), we achieve an rms of 0.73 K. 

2.3. Combination procedure 

In order to combine the data from the single dish and inter- 
ferometer, the NR045 image is converted to visibilities, com- 
bined with the CARMA visibilities in the wv-plane, and the 
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TABLE 1 

12 CO OBSERVATIONS OF NGC 6946 



TABLE 2 
Properties of NGC 6946 





NR0 45 


CARMA 


Combined 


Observing time (hours) 


24 


21 


45 


Bandwidth (MHz) 


265 


100 


100 


Velocity resolution (km s ) 


2.54 


2.54 


2.54 


Beam size (") 


19.7 




2.3 X 2.3 


RMS (K) 


0.13 


0.73 


1.9 (0.11 Jy beam"') 



new uv dataset is imaged together. We f ollow the procedure 
thoroughly described in lKoda et al.l d201 lb for the imaging of 
M5 1 , and we refer the reader to that paper for the details of the 
NR045 deconvolution, combination in Mv-space, and imaging 
process. In that paper, the relative Mv-coverage of the single 
dish and interferometer observations enabled the single dish 
visibilities to be flagged beyond 4 kA, as CARMA visibilities 
existed down to this value. We keep the NR045 visibilities 
out to 10 kA in order to ensure sufficient overlap between the 
two sets of Mv-coverage. The rms of the combined cube is 
0.11 JybeairT 1 (1.9 K) using the combined synthesized beam, 
discussed below. We maintain the instrumental velocity res- 
olution of 2.54 km s~' to optimize our ability to resolve, not 
only to detect, GMCs. The observing parameters for both the 
single dish and interferometric data sets are summarized in 
Table Q] 

2.3.1. Combined synthesized beam 

Combining interferometric and single dish information is 
non-trivial, as each is imprinted with the intrinsic beam size 
with which the observation is made. Typically, when decon- 
volving interferometric data, the process (e.g., CLEAN) re- 
moves the synthesized beam pattern at the location of each 
emission peak detected in the dirty map above some thresh- 
old and replaces it with a convolution beam, which is typi- 
cally Gaussian in shape. The convolution beam in the case 
of a combined (interferometer + single dish) data set requires 
a beam comprised of weighted components from both tele- 
scopes and is effectively a superposition of two Gaussians 
with quite different FWHM sizes. Cleaning algorithms, such 
as invert within the interferometric data reduction package 
miriad, yield an output dirty beam whose point spread func- 
tion can be fitted with one Gaussian to achieve the beam size, 
but we find that this method does not lead to flux conserva- 
tion since the two-component beam does not have a simple 
Gaussian shape. 

In order to ensure that we conserve the flux in our combined 
NR045+CARMA cube, we assume that the single dish obser- 
vations accurately measure the total flux of NGC 6946. We 
calculate the flux in the single dish dirty map and compare it 
to the flux measured in the combined cleaned map; if the flux 
is conserved, these should be equal. We choose the beam size 
which leads to flux conservation by setting equal the single 
dish flux per single dish beam (with units Jy beam s l D ), both of 
which we know, and the flux in the combined dirty map per 
combined beam (with units Jy beairT^), where we know the 
flux but not the beam size: 



flux SD (Jy) flux comhi „ ed (Jy) 



beam s D (") beam. 



combined 



CO' 



(1) 



RAy2000 


20''34 m 52.3 s 


Decj20oo 


60° 09' 14" 


Morphology" 


SABcd 


Distance* 


6.08 Mpc 


Optical velocity" 


40 km s" 1 


Major diameter" 


11.5' 



In this way, we calculate the intrinsic size of the combined 
beam (2.3" x 2.3") and use this value to clean the combined 



REFERENCES: a: NAS A/IPAC Extrag alactic Database (NED); 
b: Herrmann et aljiool 
map. We refer the reader to lKoda et al.l d201 lh for a complete 
discussion of this step in the data combination procedure. 

3. RESULTS 
3.1. Imaging 

NGC 6946 is a well-known nearby spiral galaxy. Its gen- 
eral properties are listed in Table |2l and our CO imaging is 
presented in Figures Q] and [2] The combined uv data sets 
are imaged with purely uniform weighting (robust = -5), 
where a weighting scheme is applied to the visibilities which 
is inversely proportional to the sampling density (in effect, 
weighting all cells the same at the cost of higher rms noise). 
This method takes advantage of the highest possible resolu- 
tion of the array and yields a beam size of 2.3" (calculated as 
described in Section 2.3.1). 

We overlay the high resolution CO imaging on a Hubble 
Space Telescope (HST) Ha image of the center of NGC 6946 
in Figure [3] The CO emission traces the spiral arm and bar 
structures at the center of the galaxy. To the north and south- 
east of the nucleus, the CO emission and dust lanes apparent 
in the Ha image appear to be well matched. As the spiral 
arm structure consists of adjacent stellar light and dust fea- 
tures (inferred from the Ha absorption), where the CO does 
not match the absorption precisely, its peaks tend to be found 
on the edge of the absorption feature nearer to the stellar com- 
ponent. 

NG C 6946 is known to host a molecular bar (Ball et al. 
1985), and we show the velocity field of NGC 6946 in Fig- 
ure |4] to accentuate the regions where narrow spatial regions 
span large ranges in velocity. These are the offset ridges of 
the galactic bar, which largely coincide with the edges of the 
dust lanes seen in the Ho- image. 

3.2. Defining giant molecular clouds 

Given the distance of NGC 6946 and our spatial reso- 
lution of 68 pc, we expect to be sensitive to the largest 
GMCs and giant molecular associations (GMAs). The typ- 
ical s ize of a Galactic GMC is 40 pc (Scoville & Sanders 
119871) . but measurements range from less than 2 pc to over 
100 pc dSolomon et al.lll987HScoville et al.lll987l). We utilize 
CLUMPFIND dWilliams. de Geus. & Blitd 119941) . a well- 
known algorithm designed to search image cubes in three di- 
mensional (x-y-v) space for coherent emission, in order to 
decompose our d etected CO emission into individual molec- 
ular c louds (e.g.. iKoda et al.l 12009; Eg usa. Koda. & Scovilld 
HOTTI) . 

CLUMPFIND assigns emission to individual clouds by 
contouring the emission in a data cube at multiples of the rms 
noise, identifying peaks, and following the peaks to lower in- 
tensities. We find CLUMPFIND to be the simplest algorithm 
for cloud decomposition, as it requires only two user-specified 
inputs: the contour increment and the minimum contour level 
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Fig. 1 . — CO map of NGC 6946 (left) and flux assigned to 1 34 clumps by CLU MPFIND (right). Our imaging is optimized to resolve GMCs but not to detect 
less massive GMCs. 
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Fig. 2. — CO emission from NGC 6946; every other velocity channel is shown. The contour plotted represents 2cr (0.22 Jy beam -1 ). 
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Fig. 3. — The central region of NGC 6946 is shown. CO contours are overlaid on the CO imaging (left) and Ha imaging from HST (right). The molecular gas 
coincides with the spiral arm structures, namely the apparent stellar features and dust lanes. North is up, and east is to the left. 
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Fig. 5. — Clouds detected in the central region of NGC 6946. Channels are 
separated by 5 km s , and colors are coded by clump number. The black 
circles indicate the same scale on each panel. 

(i.e., where to stop), and it assumes no physical priors in the 
clump decomposition process. We find that the most believ- 
able output is returned when using twice the rms noise of the 
data cube (0.11 Jy b eam" 1 ) for both the increment an d min- 
imum contour level; Willia ms, de Geus. & Blitz! (119941) also 
strongly recommend these values for the two parameters. We 
show an example of the clouds extracted by the algorithm in 
several subsequent velocity channels in Figure [5] 

All clumps returned by the algorithm have integrated flux 
measurements greater than 3<x (0.33 Jy beam -1 ). Addition- 
ally, the velocity profile of each clump is visually inspected 
using the clplot package within CLUMPFIND, and any 
clumps with profiles which do not appear to be single enti- 
ties (i.e., apparently blended profiles) are removed from the 
sample. 

The algorithm returns the FWHM extents of each clump 
in the x-, y-, and v-directions, calculates the total two- 
dimensional areas of each clump from the pixels assigned and 
the total detected flux within the three dimensional clump vol- 
ume, and derives clump radii assuming the calculated areas 
are circular. We account for the beam size (the spatial reso- 
lution element) by subtracting in quadrature the beam size of 
a point source with the same peak temperature as the clump, 
measured out to the radius at which the increment equals the 
minimum contour level (where T=AT), from the radii calcu- 
lated using the area method (as d escribed in Appendix A of 
I Williams, de Geus. & Blitzll 1994b . We account for the veloc- 
ity resolution, or the (non-Gaussian) channel width, in our 
velocity dispersions by approximating the instrumental dis- 
persion to be half of the channel width and subtracting this 
value in quadrature from the measured velocity dispersions. 
The equations are: 



2.355 



and 



CTy = 



_ ^ chart y 



(2) 



(3) 



iRosolowskv & Lerovl d2006) define more complex treat- 
ments of these quantities which require the extrapolations of 
clouds out to K. We find good agreement between the val- 
ues for the velocity dispersion (to within 4%) calculated using 



both methods, but the size measurements diverge for clouds 
smaller than our beam size. Since our measurements of size 
and velocity dispersion are measured out to a contour equal 
to twice the rms, not extrapolated to K, we opt to use the 
Williams method to estimate the sizes and velocity disper- 
sions of the clouds. 

Following e.g., IWilson & Scoviilel (Q990), we take the un- 
certainties in our measurements to be 25% of the spatial beam 
size (17 pc) and half of our velocity channel width (1.3 km 
s~'). These values are appropriate for the largest clouds in 
our sample, but the uncertainties in the clouds with sizes and 
velocity dispersions near our resolution limits may be larger. 
Running CLUMPFIND to a lower minimum contour level 
(1.5cr) creates more detections of faint clouds, but the result- 
ing change in Xco within the individual GMCs presented in 
this paper is well within these quoted measurement errors. 

3.2.1. Derived properties 

We use the cloud measurements to calculate cloud virial 
masses in the follo wing manner. Using the d efinition of po- 
tential energy from McKee & Zweibel ( 1992) and the kinetic 
energy in terms of the one - dimen sional velocity dispersion, 
Williams. deGeus. & Blitzl (fl994l) express the virial mass as 



M*, 



5ARcrl 
a vir G 



(4) 



where AR is the circular radius of a clump, cr v is the one- 
dimensional velocity dispersion, G is the gravitational con- 
stant, and a V i r is the "virial parameter" which describes a non- 
uniform density profile. If we parameterize the density profile 



p(n=pi\ — 



we find that integrating over dM = p(r)ydydx yields 



M vir = 



3-/3 



(5) 



(6) 



Substituting the expressions in Equations 5 and 6 into 
the gravitational energy fo r a spherical body (Equation 2.32, 
iBinnev & Trem aine 2008) yields the relation 



-a vir 



3-1 
5 -2/3' 



(7) 



for y8=l, or V!> =10/9; 
also mentioned by 



It follows that for /3=0, ar v; >=l; 
and for 0=2, a vir =5/3 ( as 
iWilliams. de Geus. & BhtzTHl^l . 

If we assume that / ?=!, as many authors do, (e.g. , 
Solomon et all 119871; IRosolows kv & Lerovl 120061; 



Bolatto et al.1 120081) and G = 1/232 (using units of km 
s _1 , parsecs, and solar mass), substituting into Equation 4 
yields 

M vir = 1040 R o*, (8) 

as shown by | Solomon et al.l (Q987) and similarly by 
IWilson & Scovilld(ll990l) . We assume that the GMCs are vit- 
alized, an assumption consistent with previous studies, and 
derive virial masses using the measurements of R and cr v (cor- 
re cted for their respectiv e resolution elements) via Equation 8. 
Solo mon et al.l (119871) describe the virial mass as 
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My, 



(9) 



where f p is called a projection factor and the size of the cloud 
(S) is related to the effective radius (derived from the circular 
area of the cloud on the sky) such that R e ff = 1.9 IS. However, 
this projection factor, indicated to equal 2.9, is not explicitly 
derived and as a result does not appear in many recent papers 
on this topic. Using the formulae shown above, it is trivial to 
derive that f p =2.9 when a v i r = 10/9. 

To calculate Xqo, in effect the mass-to-light ratio for molec- 
ular clouds, we compare the virial masses of each cloud to 
the luminosities derived from each cloud's integrated CO flux. 
We calculate cloud luminosities via 



^co 



0«0 h 



(10) 



where A is the observed wavelength (in mm), Fco is the flux 
density in Jy beam -1 km s^'arcsec 2 , and 8 a and 9t are the 
beam axes (arcsec). Finally, with M V! > in units of solar masses 
and Leo in units of K km s~'pc 2 , and including a factor of 
1.36 to account for helium, we compute Xco (the CO-to-H2 
conversion factor) as follows: 



XcofcirT 2 (K km s' 1 )' 1 ] = 4.60 x 10 19 — — . 



(11) 



^co 



4. DISCUSSION 



4.1. Structure and dynamics of CO 

From Figure Q] it is clear that most of the bright CO 
emission which is included in the resolved GMCs is co- 
incident with spiral arms, but significant emission is also 
found in the inter-arm regions. The GMCs identified by 
the CLUMPFIND algorithm tend to be the largest, brightest 
clouds: those in the spiral arms. The inter-arm emission ap- 
pears to be more extended, which is consistent with it being 
less likely to display an apparently self-gravitating velocity 
profile. 

4.2. Properties of individual GMCs 

We resolve a total of 134 clouds. After discarding poten- 
tial blends, as described in Section 3.2, 120 clouds remain in 
our GMC sample. Using the equations in Section 3.2 to ac- 
count for the instrumental resolution elements, 64 clouds have 
a real velocity dispersion. We require that clouds have a ve- 
locity dispersion of at least 2.54 km s _1 (twice our estimated 
instrumental dispersion) to be included in the sample, leaving 
30 resolved clouds. The entire sample of 134 identified clouds 
is useful to examine the overall distribution of molecular gas 
which exists in GMCs, but the 30 fully resolved clouds will 
be the ones for which we derive virial masses, luminosities, 
and Xqo- These clouds are listed in Table [3] 

The radii and velocity dispersions of the 30 fully resolved 
GMCs are shown in Figure [6] as are the m easurements of 
resolv ed Milky Way disk clouds detec t ed by ISolomon et al.l 
(fl987h for comparison. Solomon et al. (1987) performed the 
seminal study of GMCs in the Galactic disk and showed that 
the relationship between the sizes and velocity dispersions of 
GMCs are not related by an exponent of |, as predicted by 
Kolmogorov turbulence dLarsonlfl981h . but instead by an in- 
dex of 0.5, as is consistent with clouds in virial equilibrium. 



However, more recently, comparisons of GMCs at the 
Galactic center using various tracers have been made to those 
throughout the disk, indicating that departures from the tra- 
ditional relation do occur for clouds at the Galactic center 
dMiyazaki & Tsuboil 120051 lOka et all 1200 lh . The observed 
tendency is for clouds to exhibit higher velocity dispersions 
for a given size. Although our map covers the central re- 
gion of NGC 6946, the measurements of its GMCs are largely 
consistent with the Galactic disk relation with some scatter, 
particularly for the largest clouds. The measurements of ra- 
dius and velocity dispersion may even be consistent with a 
slope larger than 0.5, the value measured by several authors 
for Milky Way clouds w i thin the disk and at the Galactic cen- 
ter (ISolomon et al. 111981: fMivazaki & Tsuboil2000HOka et alj 
2001). 

For the largest clouds, the observed scatter goes to higher 
cr,, for a given size, and in fact, the six clouds which fall well 
above the Solomon relation are all located within ~400 pc of 
the center of NGC 6946, consistent with the trend observed 
to be true for the center of our Galaxy. The smallest clouds 
which deviate from the Solomon relation are very close to our 
velocity resolution limit. While we are not sensitive to the 
full ra nge of GMCs detected in the Galaxy by Solomo n" et alJ 
(1987), our detected GMCs are largely consistent with the 
biggest clouds in the Galactic sample. Whi le the Milky Way 
clouds included in the ISolomon et al.l ( 1 19871) study were mea- 
sured using a contouring method which is physically distinct 
from the decomposition algorithm employed in this paper, the 
consistent trends present in both samples indicate similar un- 
derlying physics at work in both the center of NGC 6946 and 
the disk of our Galaxy. 

The GMC luminositie s are shown in Figure|7] The Galactic 
points and their linear fit ( Solomon et al.ll987l) are also shown 
for comparison. Again, our GMC sample is largely consistent 
with the measurements of the brightest Galactic clouds in the 
Solomon sample, though a fit to the NGC 694 6 clouds would 
be sha llower than the relation measured by Solomon et al. 
(1987). The brightest GMCs - the same clouds which fall 
above the Solomon radius-cr,, relation - appear to be underlu- 
minous for their velocity dispersion compared to the Galactic 
relation. Though blended clouds could artificially produce 
this relationship, since the total luminosity of multiple op- 
tically thick 12 CO clouds along the same line of sight may 
not increase linearly with the amount of CO-emitting gas (de- 
pending on the geometry), the velocity dispersion would be- 
tray the presence of more than one cloud via multiple peaks. 
We individually inspect each cloud to select them on the ba- 
sis of their apparently unblended velocity profiles, so blended 
clouds (at least at our velocity resolution limit) can not explain 
this finding. 

4.3. CO-to-H^ conversion factor 

The virial masses and CO luminosities of our NGC 6946 
GMC sample are shown in Figure[8] Once again, we compare 
our measureme nts to the Galac t ic sa mple; the overlaid linear 
fit is derived in ISolomon et alj d 1987b . The median mass of 
our GMC sample is comparable to the most massive Galac- 
tic disk GMCs, but we find excellent consistency between the 
trends. The histogram of our GMC sample is shown in Fig- 
ure!] 

Including a factor of 1.36 to account for heavier elements, 
the average value of Xco for the GMC sample is 1.2 x 10 20 
cirT 2 (K km s _1 ) ; this value can vary by as much as a factor 
of 2 with our assumed uncertainties. This value is roughly a 
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tic Solomon et al. 1 1987) samples are shown. The fit to the Galactic points, 
L co ~ 130 o-^ Kkms" 1 pc 2 , is also shown. 



factor of two below the averag e value of the conve rsion factor 
which we calculate from the Solomon et al. (1987) Galactic 
disk clouds of2.6x 10 20 cirT 2 (Kkms -1 ) -1 (derived using the 
sa me factor for heavier e lements). If we calculate Xco from 
the lSolomon et al.l (|T987) clouds using only GMCs more mas- 
sive than 5 x 10 5 M Q - approximately our completeness limit 
- we find a value of 1.9 x 10 20 cm (K km s -1 ) -1 , which is 
in agreement with the value we derive for NGC 6946 within 
our unc ertainty. It is also in agreement with the value de- 
rived bv lDame, Hartmann, & Thaddeusl ( 1200 11) of 1.8 (+0.3) 
x 1 20 cm -2 (K km s~ 1 )~ 1 for CO emission in the solar neigh- 
borhood. 

In addition to CO measurements, the conversion factor in 
the Galaxy has also bee n derived using y-ray measurements 
dStron g & Mattoxl [l996l) . which yield a conversion factor of 
1.9 x 10 20 cirT 2 (K km s -1 ) -1 , in good agreement with the 
values listed above. This agreement between values of Xco 
determined using both dynamical and non-dynamical meth- 
ods has been claimed as evidence that Galactic GMCs are 
gravitation ally bound and in eq uilibrium. iHever et ail (120091) 
correct the Sol omon et al.l (11987b measurements for a more re- 
cently derived Galactic center distance and a factor of 1.36 to 
account for heavier elements and derive a mean surface den- 
sity of 206 M Q pc~ 2 for the entire Galactic disk GMC sample; 
the mean surface density that we calculate within our clouds 
is 204 M pc~ 2 . 

The scatter in measurements of the conversion factor among 



studies of nearby galaxies has tended historically to be large, 
but more recently, studies of Xco ha ve exhibited relativ e 
agreement with the Galactic value. In Bol atto et al.l (2008), 
12 CO (J=l-0) and (J=2-l) measurements of GMCs within Lo- 
cal Group dwarf galaxies and spiral galaxies are compiled and 
compared to the molecular clouds detected in the Milky Way. 
The sample includes dwarf galaxies as well as the Large and 
Small Magellanic Clouds (LMC and SMC) to investigate the 
effect played by the environment on determinations of .Xco- 
These authors find that GMCs in Local Group dwarfs obey 
similar relationships to Milky Way GMCs and are consistent, 
to within a factor of two, with an Xco of 2 x 10 20 crrT 2 (K km 
s -1 ) -1 with no detectable dependence on metallicity. A com- 
parison of the GMC masses and luminosities of the NGC 6946 
sampl e and the extragalactic G MC sample f rom Bolatt o~et alj 
(2008) is presented in FigureflOl B litz et all (120071) also study 
GMCs within nearby galaxies and find that Xco is within a 
factor of two of 4 x 10 20 cirT 2 (K km s -1 ) -1 . Even GMCs in 
the outer disk of M33 exhibit charact eristics similar to those 
of Galactic GMCs dBigiel et aUl2010t). I n both the studies by 
iBolatto et all d2008l) and bv lBlitz et all (120071). howev er, the 
low-metallicity SMC is an outlier: in lBlitz et al.l (120071) . it re- 
quires a conversion factor three times higher than the rest of 
the GMCs in the sample. 

As it is common for disk galaxies to exhibit a metallic- 
ity gradient, we investigate whether a radial dependence of 
the distributions of sizes, masses, and Xco values of individ- 
ual GMCs exists. Various observational studies which cor- 
relate the conversion factor with metallicity do so by study- 
ing them as a function of radius (or even averaged over en- 
tire disks) to investigate how Xco may vary with metallicity 
in the presence of a metallici ty gradient, especially in dif- 
fuse and low-met allicity gas (Arimoto . Sofue. & Tsuiimotol 
| 1996t iLerovet al.l 1201 lb iGenzel et al.l 120111) . For instance. 
iNakai & Ku no (1995) show that M51 exhibits a radial gradi- 
ent in the conversion factor in M5 1 , with Xco increasing with 
radius as the metallicity decre ases. The same has also been 
shown to be true in the Galaxy (Sodroski 199Tb. 

These properties are displayed in Figure [TT] In this figure, 
the full detected cloud sample of 134 clouds is shown as dots. 
The 120 clouds with velocity profiles which point unambigu- 
ously to self-gravitation are circled in a color which indicates 
the mass (left) or value of Xco (right) while the size of the 
circles indicates the radius of each cloud. The 30 clouds used 
in our analysis are indicated by asterisks. The most massive 
clouds (> 10 7 M ) are found at the very center of NGC 6946, 
as described in Section 4.2, and clouds with masses > 5 x 10 5 
M are found where the density of clouds is the highest (i.e., 
on the spiral arms). No radial trend is apparent in Xco, though 
the most central clouds tend to have values higher than 0.5 x 
10 20 cirT 2 (K km s _1 ) , and all but one of the clouds with 
Xco below this value are smaller than 70 pc. 

4.3.1. Metallicity 

At radii less than 2.5 kpc from the center of NGC 6946, the 
mean O/H abundance level is roughly constant ( B ellev & Rovl 
119921) ; this is the region over which we present resolved mea- 
surements of GMCs, so we do not expect systematic or radial 
changes in metallicity to affect our measured values of .Xco- 
However, it is still worthwhile to compare our derived value 
of the conversion factor to the metallicity at the center of NGC 
6946 in the context of previously derived trends relating the 
two. Typically, the conversion factor is found to increase as 
the host galaxy metallicity decreases. The metallicity mea- 
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Fig. 8. — Lum inosities and virial m asses of the NGC 6946 GMC sample 
and the Galactic Solomon et al. ( 1987) sample are shown. The overplotted fit 
to the Galactic points. M,,„ = 39 X (Leo) 81 , is consistent with both sets of 
GMCs. 

sured at the center of NGC 6946 quoted by Moustaka s et al.l 
(2010) is [12 + log (O/H)] = 8.47+0.09 or [12 + log (O/H)] 
= 9.16+0.06, depending on the empirical or theoretical abun- 
dance calibration applied. The two values are intended to in- 
dicate the range of metallicities which would be derived using 
current strong-line calibrations (Moustakas et al. 2010). 

Our derived value of Xco and the lower Moustakas value of 
the metallicity are consistent with the relation a/acai = 6.0 - 
0.67[12 + log (O/H)] found bv lWilsori ( fl995l) . where the cor- 
relation between Xco and metallicity is derived using virial 
GMC measurements. In addition to using virial measure- 
ments, the metallicity dependence of Xc o has also been teste d 
by assuming that dust traces gas (i.e., Ma grini et al.M2.0lTb . 
Our derived value of Xco and the higher Moustakas value 
of the metallicit y are consi sten t with the trends found using 
this m ethod in llsraell (120001) and lBoselli. Lequeux. & Gavazzil 
(2002), which are s hown over 2-3 orde rs of magnitude in .Xco- 
On the other hand, iLerov et ail (1201 ll) find that metallicity is 
correlated with the conversion factor at metallicities below 
[12 + log (O/H)] = 8.2; the entire Moustakas range of val- 
ues for the metallicity at the center of NGC 6946 is higher 
than this value. 

In the Milky Way, iBalser et all d20lTI) measure a radial 
metallicity gradient in a sample of Galactic disk HII regions 
to be [12 + log (O/H)] = (8.962 + 0.045) - (0.0446 + 0.0046) 
R ga i, which yields Galactic disk metallicities in the range of 
8.5 to 8.7 at radii of 5-10 kpc. These measurements of the 
Galactic disk metallicity are consistent with the Moustakas 
range of values of the metallicity at the center of NGC 6946; 
thus, it is not unreasonable to expect that the GMCs in the two 
regions are physically similar, especially within our quoted 
uncertainty for Xco- 

4.4. "Diffuse " fraction of CO 

The fraction of emission in our CO cube which is identified 
as emission from the 30 GMCs with approximately Gaussian 
velocity profiles (suggesting self-gravitation) greater than our 
velocity resolution element is 19%; the total mass in this com- 
ponent is 2.2xl0 8 M Q . If we include all of the GMCs identi- 
fied by CLUMPFIND, which includes clouds which are spa- 
tially resolved but exhibit blended velocity profiles, the frac- 
tion increases to 31%. However, since the algorithm discards 
emission contours below 2<x, fainter emission around GMCs 
is not included. While this does not affect Xco, as described 
in Section 3.2, it does artificially elevate the "diffuse" fraction 
that we find. For instance, the fraction of emission included 



Fig. 9. — Histogram of the GMC masses in the NGC 6946 sample. The 
solid line indicates the sample of 120 resolved clouds with apparently self- 
gravitating profiles, and the shaded area indicates the sample of 30 clouds on 
which we base our analysis. 
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Fig. 10 . — The same as Fi gure [8] Extragalactic GMC properties from the 
study by Bolatto et al. (2008), and their corresponding fit, are shown as well 
as the[Solomon et al. 1 1987) fit to the Milky Way sample. The clouds mea- 
sured in NGC 6946 are more consistent with the fit to the Galactic GMCs. 

in all GMCs in the 1.5<xrun is 47% (compared to 31%). 

The critical density for collisional excitation of CO 
(J=l-0) emission is 300 enr 3 dScoville & Sanders! 11987b . 
which is roughly the average gas density within a GMC 
(ISolomon et al.l[l987l) . The CO emission that is not identi- 
fied as belonging to GMCs in our analysis is unlikely to be 
entirely diffuse (i.e., below the critical density). In addition, 
our velocity resolution of 2.54 km s _1 - optimal for resolving 
more massive GMCs - is not optimal for the detection of less 
massive GMCs. 

We expect to be sensitive to clouds down to a limit of 
~5 x 10 5 M , but certainly clouds less massive than this 
lim it will contribute as well. In their Galactic GMC sam- 
ple, [ScovmiAlindiri 0113) find that -50% of the total H 2 
mass resides in clouds less massive than 4 x 10 s M Q , and the 
other half is found in clouds more massive than this value. 
Taking our estimate - that roughly half of the CO emission 
that we detect (in our 1 .5cr run of CLUMPFIND) is assigned 
to cloud s more massive than ~5 x 10 5 M Q - in conjunction 
with the lScoville & Sandersl (119871) result for the Galaxy in- 
dicates that the fraction of truly diffuse CO that we detect is 
likely small. 

5. CONCLUSIONS 

In this paper, we present the first galaxy from our CO Sur- 
vey of Nearby Galaxies taken with CARMA and NR045. 
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Fig. 1 1 . — Spatial distribution of radii, masses, and Xco values. In these plots, all 1 34 resolved clouds are shown as dots. The 1 20 clouds whose velocity profiles 
indicate that the clouds are apparently self-gravitating entities are identified by a colored circle as well. The 30 clouds which have real velocity dispersions after 
subtracting the velocity resolution element, which used in our analysis, are indicated by asterisks. The radii, mass, and Xco of each cloud with such a velocity 
profile is coded by color and size to indicate the values of these properties. 



The high image fidelity afforded by our combined observa- 
tions allows us to study extragalactic molecular clouds with 
unprecedented resolut ion. We utilize the clump-fi nding algo- 
rithm CLUMPFIND dWilliams. de Geus. & Blitzll 19941) with 
the recommended rms noise-dependent inputs and find that 
the general tren ds recovered by the cl ouds resemble those of 
Galactic clouds dSolomon et al.|[l987l) . but the clouds within 
~400 pc of the center of NGC 6946 exhibit higher velocity 
dis persion fo r a giv en size, similar to the behavior observed 
by lOka et al.l ((2001) for clouds at the Galactic center. These 
clouds also exhibit higher velocity dispersions for a given lu- 
minosity, which is not likely to be caused by blends of clouds 
(at least at our velocity resolution limit). 

We derive an average value of Xco of 1.2 x 10 20 crrT 2 
(K km s -1 ) -1 , which can vary by as much as a factor of 2 
with our assumed measurement uncertainties, and is consis- 
tent within our errors with accepted Galactic values of Xco 
(Dame, Hartma nn. & Tha ddeus 2001). The trends which we 
observe among the GMCs at the center of NGC 6946 are 
broadl y consistent with Galactic trends seen by Solomon et al. 
(119871) . even though we use a phenomenologically different 
algorithm to define the clouds, and we detect clouds compa- 
rable to the brightest and most massive clouds in the Galactic 
sample. The most massive clouds (> 10 7 M ) are found at the 
very center of NGC 6946, and clouds with masses > 5 x 10 5 
M are found where the density of clouds is the highest (i.e., 
on the spiral arms). No radial trend is apparent in Xco, though 
the most central clouds tend to have values higher than 0.5 x 



10 20 crrT 2 (K km s -1 ) -1 , and all but one of the clouds with 
Xco below this value are smaller than 70 pc. 

When we run the clump decomposition algorithm down to 
1.5cr, which does not change our measured value of Xco but 
does incorporate more extended emission around our detected 
GMCs, we find that 47% of our detected CO emission is iden- 
tified as belonging to a GMC. Our study is optimized to re- 
solve larger and more massive GMCs as opposed to detecting 
less massive GMCs; we ar e most sensitive to clouds a bove ~5 
X 10 5 M . If we adopt the lScoviUe & Sanders! (|1987l) finding 
in the Galaxy that 50% of the H2 mass is in clouds more mas- 
sive than 4 x 10 5 M and 50% is in clouds less massive than 
this value, then the fraction of truly diffuse CO emission that 
we detect in NGC 6946 must be small. 

Studies of this nature will soon be expanded when ALMA 
comes online, as observations will be feasible with better sen- 
sitivity and resolution in a fraction of the observing time. 

The authors would like to thank the anonymous referee for 
the suggestions which improved the discussion of this paper. 
This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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TABLE 3 
Properties of GMCs in NGC 6946 



Number 


RA/2000 


Decj2000 


R(pc) 


cr v (km s 1 ) 


L co (10 5 Kkms-'pc 2 ) 


M vir (10 s Mo) 


Xqo (10 20 [cm" 2 (Kkm s" 1 )-']) 


1 


20* 34 m 52.73 s 


60° 9' 12" 


117 


11.5 


118 


161 


0.630 


2 


20* 34 m 52.66 s 


60° 9' 14" 


155 


14.7 


152 


348 


1.05 


3 


20* 34™ 52.73 s 


60° 9' 13" 


154 


7.89 


88.5 


99.7 


0.518 


4 


20* 34 m 53.60 s 


60° 9' 19" 


189 


9.04 


117 


161 


0.633 


5 


20* 34 m 51.79 s 


60° 9' 18" 


123 


8.86 


55.8 


101 


0.829 


6 


20* 34 m 52.73 s 


60° 9' 


89.7 


10.5 


31.7 


104 


1.50 


7 


20* 34 m 52.80 s 


60° 9' 30" 


146 


9.91 


46.2 


149 


1.49 


8 


20* 34™ 48.05 s 


60° 7' 53" 


145 


5.84 


28.9 


51.6 


0.822 


9 


20* 34 m 51.66 s 


60° 9' 23" 


113 


14.2 


50.1 


236 


2.17 


10 


20* 34™ 41.48 s 


60° 8' 47" 


43.8 


3.82 


3.09 


6.64 


0.989 


11 


20* 34 m 52.20 s 


60° 9' 27" 


94.8 


12.0 


28.7 


143 


2.29 


12 


20* 34 m 51.93 s 


60° 9' 19" 


95.8 


17.6 


42.0 


309 


3.38 


13 


20* 34 m 45.77 s 


60° 8' 


138 


5.52 


24.4 


43.7 


0.825 


14 


20* 34 m 41.21 s 


60° 8' 49" 


40.4 


4.89 


4.59 


10.0 


1.01 


15 


20* 34 m 42.95 s 


60° 9' 6.5" 


116 


8.01 


23.0 


77.1 


1.54 


16 


20* 34 m 50.72 s 


60° 8' 25" 


98.0 


5.87 


13.5 


35.1 


1.20 


17 


20* 34 m 55.08 s 


60° 9' 51" 


45.2 


4.90 


4.66 


11.3 


1.11 


18 


20* 34 m 43.22 s 


60° 8' 35" 


76.0 


2.76 


6.49 


6.04 


0.428 


19 


20* 34 m 59.76 s 


60° 8' 18" 


89.5 


5.02 


8.55 


23.4 


1.26 


20 


20* 35™ 2.913 s 


60° 9' 13" 


57.8 


2.63 


3.33 


4.17 


0.576 


21 


20* 34 m 55.88 s 


60° 8' 25" 


60.3 


2.76 


3.22 


4.79 


0.683 


22 


20* 34™ 56.01 s 


60° 8' 15" 


82.0 


3.84 


5.87 


12.6 


0.986 


23 


20* 34™ 58.09 s 


60° 8' 37" 


37.4 


2.56 


2.06 


2.55 


0.569 


24 


20* 34 m 53.80 s 


60° 8' 10" 


49.3 


2.85 


2.21 


4.15 


0.863 


25 


20* 34 m 55.34 s 


60° 8' 53" 


59.3 


4.59 


4.06 


13.0 


1.48 


26 


20* 34 m 54.67 s 


60° 8' 5.0" 


56.9 


5.02 


3.75 


14.9 


1.83 


27 


20* 34 m 50.92 s 


60° 9' 58" 


90.3 


4.19 


6.66 


16.5 


1.14 


28 


20* 34 m 53.13 s 


60° 10' 


65.7 


3.51 


3.80 


8.42 


1.02 


29 


20* 34 m 50.32 s 


60° 8' 27" 


60.9 


3.15 


3.43 


6.29 


0.844 


30 


20* 34 m 53.07 s 


60° 8' 59" 


39.2 


5.40 


1.89 


11.9 


2.90 
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